Landauer-type transport theory for interacting quantum wires: Application to carbon 

nanotube Y junctions 
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We propose a Landauer-like theory for nonlinear transport 
in networks of one-dimensional interacting quantum wires 
(Luttinger liquids). A concrete example of current experi- 
mental focus is given by carbon nanotube Y junctions. Our 
theory has three basic ingredients that allow to explicitly solve 
this transport problem: (i) radiative boundary conditions to 
describe the coupling to external leads, (ii) the Kirchhoff node 
rule describing charge conservation, and (iii) density matching 
conditions at every node. 

PACS numbers: 71.10.Pm, 72.10.-d, 73.63.-b 



The Landauer-Biittiker approach to transport in meso- 
scopic systems has been very successful in describing non- 
interacting electrons by using a scattering matrix for- 
mulation ||l[. It is an important challenge to general- 
ize this approach to the case of strongly correlated elec- 
trons. Here we propose such a theory for iV-terminal star- 
like networks of interacting ID quantum wires (QWs) 
described by Luttinger liquid (LL) theory. The two- 
terminal setup (A^ = 2) has been formulated and solved 
previously [^js). However, the step from A^ = 2 to A^ = 3, 
where three individually contacted QWs meet at a single 
node ("Y junction"), see Fig.|^, is non-trivial yet essen- 
tial for the development of a practically useful transport 
theory for interacting electrons. Recently, several works 
appeared where precisely this problem has been under 
study. While in Ref. |^ a weakly coupled "Kondo" node 
was considered, other authors used perturbation theory 
in the hopping j|] and/or the interaction At the same 
time, it has become clear that a more general approach 
is necessary to go beyond those special situations. Below 
we formulate boundary conditions that allow for the ex- 
plicit solution of this transport problem. Progress along 
these lines is also likely to sharpen our understanding of 
conformal field theory with boundary conditions . 

This problem is not only of intellectual interest but 
also of relevance to transport experiments for carbon nan- 
otubes [^. As has been predicted j|] and observed in a 
series of beautiful experiments , single- wall nanotubes 
provide a realization of LL physics. Electron-electron in- 
teractions cause remarkably pronounced non-Fermi liq- 
uid behaviors characterized by the standard LL parame- 
ter g « 0.25 (where g = 1 is the non-interacting value). 
Template-based chemical vapor deposition |ll] and elec- 
tron beam welding methods [|2| allow to fabricate and 



contact nanotube Y junctions, and the intrinsic nonlinear 
I — V characteristics of such a device has been observed 
recently jl3j. In addition, the case A^ = 4 has been real- 
ized by several groups using two crossed nanotubes [ p^ , 
providing another interesting application. Furthermore, 
semiconductor heterostructures |15| or ultracold trapped 
atomic gases |l6|l may allow for the systematic study of T- 
or Y-junctions as well. Eventually it could even be possi- 
ble to access the fractional statistics of LL quasiparticles 
through the noise properties of such a device, thereby re- 
alizing a Hanbury-Brown-Twiss setup for fractionalized 
quasiparticles 

We study A^ single-channel spinless QWs at —L < x < 
described by LL theory merging at a common node 
at a; = 0. For simplicity, we assume the same interac- 
tion constant g and Fermi velocity in each QW, with 
straightforward generalization also to include spin and 
flavor degrees of freedom [|| . Theory then has to address 
(i) the inclusion of applied voltage sources and (ii) how a 
consistent treatment of the node can be achieved. Let us 
start with the first issue. Like in the two-terminal case, 
adiabatically coupled external reservoirs held at electro- 
chemical potentials fit lead to radiative boundary condi- 
tions 1^ (we put h — I) 
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where Pi^R/L{x) is the right/left-moving part of the den- 
sity in QW i. These boundary conditions only depend 
on the current injected into each QW from the respective 
reservoir, which in turn does not depend on the backscat- 
tering happening later on at the node or within the QW. 
Here we discuss the computation of the nonlinear con- 
ductance matrix normalized to e^/Zi, 
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where the current Ii flowing through QW i is oriented 
towards the node, leaving noise properties to a future 
publication. Below applied voltages Ui are deflned as 
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Under this deflnition, gauge invariance (conductance ma- 
trix is invariant under a uniform change of all ^i) is au- 
tomatically fulflUed if the depend only on the Ui . 
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Next we address the physics arising at the node x = 0. 
Conservation of charge enforces the Kirchhoff node rule 



(4) 



A second requirement at the node is the wavefunction 
matching via the 5* matrix ||l7| , |l8|l 



(5) 



where 'i'{x) = (jpi, . . . ,^Pn) contains the wavefunctions 
for the N QWs, and the outgoing (left-moving) compo- 
nents L are connected to the incoming (right-moving) 
states ^'i? via an appropriate N x N matrix S, see, 
e.g. Ref. for specific choices at TV = 3. Note that 
the scattering matrix in Eq. (||) provides a "bare" (micro- 
scopic) description of the node properties, while interac- 
tions could dynamically generate some different bound- 
ary condition at low energy scales. 

Unfortunately, a boundary condition like Eq. (|^) is 
very difficult to handle for correlated electronic systems 
and typically does not allow for progress. Here we pro- 
ceed differently by constructing a wide class of practi- 
cally important S matrices (albeit not all possible ones 
1^ ) in the following way. We first consider an ideal sys- 
tem composed of impurity-free QWs symmetrically con- 
nected at the node. Microscopically, the corresponding 
non-interacting problem could be modelled as N tight- 
binding chains with equal hopping matrix element to, 
where the "last" site of each chain is connected to the 
common node site via the same to and all on-site ener- 
gies are equal. Such a node corresponds to the special 
highly symmetric S matrix 
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(2 - N)/N J 



For this S matrix, Eq. (^) directly implies ipi{0) = ■ ■ ■ = 
iPn{0) for the components of 'J = ^'^ + Since phases 
in V'i(O) can be gauged away, a density matching condi- 
tion upon approaching the node results, 



(7) 



where pi{0) denotes the total electronic density in QW 
i close to the node. Remarkably, the conditions (|l|), (^) 
and (^ then allow to explicitly solve this transport prob- 
lem for arbitrary interactions because the condition (|^) 
does not involve wavefunction phases but only ampli- 
tudes. To arrive at more general S matrices, in a second 
step we then consider additional impurities in the differ- 
ent QWs displaced slightly away from the node. Inclusion 
of impurities does not cause conceptual difficulties, and 
such a modelling allows to construct all S matrices of 
practical interest. 



Let us first discuss the ideal node defined by the bare 
scattering matrix (||). For g = 1, the conditions (|l|), (^) 
and (Q) are in fact equivalent to the standard Landauer- 
Biittiker approach. This is easily seen by using the usual 
scattering states, e.g. the state injected into QW i — 1: 
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with reflection (transmission) amplitude r (t). For 5 = 1, 
the boundary conditions (|^) are equivalent to a Fermi 
distribution for occupying the states (|^), while Eq. (^) 
gives 1 — r + {N — l)t. Furthermore, the density match- 
ing condition (Q) yields |l-|-r|^ = \t\^ . Combining both 
equations immediately gives r — {2 — N)/N and t — 2/N, 
and hence reproduce Eq. (^. The conductance matrix is 
then 



l-il-2/Nf, G,^, = ~{2/Nf 



(9) 



Will the conductance matrix for this ideal case be af- 
fected by interactions (17 < 1)? Based on the two- 
terminal case [||, one might suspect that there is no 
renormalization because of the Fermi-liquid character of 
the leads, and hence Eq. ^ would stay valid for arbi- 
trary g. However, it turns out that for small applied 
voltages and low temperatures, the system always flows 
to the disconnected-node Rxed point, 



= 0. 
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Thus the only stable generic fixed point of this system for 
any g < 1 represents an isolated node weakly connected 
to N broken-up QWs, even for an arbitrary S matrix of 
the node. The corrections to Eq. ( p^ due to finite Ui or 
T are then sensitive to interactions. This phenomenon is 
a consequence of the strong correlations in the LL, which 
here induce asymptotically vanishing currents even in a 
perfectly clean (impurity-free) system. Equation ( [l0| ) 
also implies that open boundary bosonization |po|| allows 
to access the asymptotic low-energy transport properties 
of a QW network. 

For arbitrary g and sufficiently far away from the node 
such that Friedel oscillations ||2l| have decayed and the 
boundary conditions (|l|) hold, the left- and right-moving 
densities must combine to give 




Pi = Pi,R + Pi,L 



Here the matrix has been defined whose entries de- 
pend on g and all Ui. The reduce to standard 
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Landauer-Biittiker transmission probabilities for g = 1, 
but in general cannot be interpreted as single-particle 
quantities. It is important to stress that Eqs. (|ll|) and 
( p^ ) are consistent with both the LL equation of motion 
and the boundary conditions (|l|) for arbitrary Tij . The 
Kirchhoff node rule can then be satisfied by requiring 
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and we now use the density matching conditions (|^) to 
obtain the Tij and hence the conductance matrix. We 
mention in passing that the usual relation Tij = 1 
should not be enforced since gauge invariance under uni- 
form potential shifts in all reservoirs has been encoded in 
Eq. already. 

At this point it is crucial to realize that close to the 
node, the density Pi{x) will deviate from the spatially 
homogeneous value pt in Eq. due to Friedel oscilla- 
tions. This happens already for g = 1, as can be easily 
checked from Eq. (||) and the subsequent solution. The 
total density very close to the node is Pi{0) = pi + Spi{0), 
where 5pi{0) is the Friedel oscillation amplitude at the 
node location. The 2kF oscillatory Friedel contribution 
Spi{x) in QW i arises due to interference of the incoming 
right-movers and the left-movers that are backscattered 
at the node. Importantly, left-movers that are transmit- 
ted from the other — 1 Q Ws into QW i cannot interfere 
with the incoming right-mover and will therefore not con- 
tribute to dpi{x). This implies that 5pi{x) is identical to 
the corresponding Friedel oscillation in a two-terminal 
setup with the same (bare) reflection coefhcient as the 
one induced by the clean node, R ~ {\ — 2/iV)^. Fortu- 
nately, this allows us to obtain 5pi{Q) in an exact manner 
using the relation [|j 
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where the "four-terminal voltage" parameter Vi is found 
from a self-consistency equation. This equation is explic- 
itly given and solved for arbitrary g in Ref . |^ , and using 
this solution, we also have the pinning amplitude ( p^ ) of 
the Friedel oscillation for any value of g. The relation to 
the two-terminal problem is not a necessary ingredient 
for our calculation, but quite convenient in allowing to 
compute Spi{0) from known results. 

For concreteness, in what follows we consider g = 1/2 
where this self-consistency equation is quite simple {ks = 
1), 
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where ^p is the digamma function and Tb = ttX^/uJc for 
bandwidth lUc and impurity strength A of the correspond- 
ing two-terminal problem. (Specifically, for = 3, to 



match the reflection coefficient of the ideal Y junction, 
(ttX/uJc)'^ = 1/(1 + N/{N -2)). The solution to Eq. (|l|) 
then yields Vi as a function of Ui alone, which, however, 
itself depends on all the chemical potentials, see Eq. (|^). 
The density matching conditions ^ are then solved by 
enforcing 

a + Tn)U,-2V,=n,U, 

for all pairs k ^ j ^ 1, ■ ■ ■ N , and with the Kirchhoff rule 
we find 
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Note that the depend only on the applied voltages Ui 
but not on the pi. As a result, gauge invariance is en- 
sured. The relations (|l^) and ( p^ represent the complete 
solution for the special S matrix (^. 

Before proceeding to more general S matrices, let us 
discuss these results for the ideal Y junction (A^ ~ 3) at 
g =1/2. From Eq. (||), we get 
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where i ^ j ^ k 'm the second equation. Note that the 
conductance matrix is symmetric. For eUi <^ T, the 
linear conductances follow: 



(2S,j - 2/3) 



1 - cip'{c+ 1/2) 
l + ci^'{c+l/2)' 



where c = Tb /2ttT and i/j' is the trigamma function. As 
T — > 0, the conductance matrix approaches the stable 
fixed point (pO|), dj ~ (T/TbY- In general, for 5 < 1, we 
find Gy - {T/TbY'^^'^ as T 0. Since this power law 
coincides with the prediction of open boundary bosoniza- 
tion around ([lo|), this also provides a consistency check 
for our calculation. Corresponding power laws also gov- 
ern the nonlinear conductances. This is shown in Fig. |^ 
for g = 1/2, depicting Gii(C/) at different temperatures. 
In the T = limit, Gn ~ U^. Clearly, interactions have 
a rather spectacular effect on the transport properties of 
this system. 

Next we briefiy outline how to construct more general 
node iS* matrices based on this solution of the ideal junc- 
tion for arbitrary g. The idea is to add impurities of 
strength in each QW close to the node, Xi « —1/kp, 
which will modify the bare S matrix; for explicit calcu- 
lations, see Ref. |l^. One can then compute the G^ 
for this more general case, but still allowing for arbitrary 
g < 1, e.g. by using perturbation theory around the above 
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solution of the ideal junction. Focusing on TV = 3 and 
just one impurity, A2 = A3 = 0, to lowest order in Ai, a 
straightforward calculation gives li = /" + 5Ii , where /f 
is the current through QW i for A^ = discussed above, 
and 

A^ 

5Ii = -2eg+— sin(7rg+) cos{-Kg+) (18) 

UJc 

X r(l-25+)(2^/i7ec^e)29+-i, 

where g+ — '^g/i and 5Ik^i = —5Ii/2. Obviously, this 
perturbative estimate breaks down at very small energy 
scales for g+ < 1 but is valid for all energies at g+ > 1. 
It is straightforward to perform similar calculations for 
more than one impurity, other N , and/or higher orders in 
the Aj. From Eq. ( p^ ) and generalizations, we infer that 
Gij = is indeed the only stable fixed point. It is also 
clear that for at least one very strong impurity Ai, the 
system will reduce to one of the special cases considered 
previously [^HJ^. For g close to 1, we can also make 
explicit contact to Ref. |^ . For S matrices close to Eq. (^) 
and weak interactions, our solution indeed reproduces 
the results of Rcf . |6j . Note that the restriction to weak 
interactions allows to easily treat more general S matrices 

To conclude, we have proposed a Landauer-type theory 
for strongly interacting electrons in branched quantum 
wires such as carbon nanotube Y junctions. A broad class 
of S matrices can be covered by formulating a suitable 
boundary condition ( "density matching" ) to describe an 
ideal symmetric junction, and on top adding effective im- 
purities in the individual wires. The only stable fixed 
point corresponds to disconnected quantum wires, with 
corrections revealing the Luttinger liquid physics via var- 
ious power laws. 
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FIG. 1. Schematic view of a Y junction. QWs extend from 
adiabatically contacted reservoirs with electro-chemical po- 
tential Hi at X — to the node at a; = 0. 




FIG. 2. Nonlinear conductance Gii((7) for A'^ = 3, 
fii = Ef + U , IJ.2 ~ I-I3 ~ Ef, several temperatures and 

g = 1/2. 
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